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ABSTRACT 

We present 21-cm observations and models of the H 1 kinematics and distribution of NGC 4244, a 
nearby edge-on Scd galaxy observed as part of the Westerbork Hydrogen Accretion in LOcal GAlax- 
ieS (HALOGAS) survey. Our models give insight into the H 1 kinematics and distribution with an 
emphasis on the potential existence of extra-planar gas as well as a negative gradient in rotational 
velocity with height above the plane of the disk (a lag). Our models yield strong evidence against 
a significantly extended halo and instead favor a warp component along the line of sight as an ex- 
planation for most of the observed thickening of the disk. Based on these models, we detect a lag of 
—9^2 km s _1 kpc -1 in the approaching half and — 9±2 km s _1 kpc^ 1 in the receding half. This lag 
decreases in magnitude to — 5±2 km s _1 kpc -1 and — 4±2 km s" 1 kpc" 1 near a radius of 10 kpc in the 
approaching and receding halves respectively. Additionally, we detect several distinct morphological 
and kinematic features including a shell that is probably driven by star formation within the disk. 
Subject headings: galaxies: spiral - galaxies: halos - galaxies: kinematics and dynamics - galaxies: 
structure - galaxies: individual (NGC 4244) - galaxies: ISM 



1. INTRODUCTION 

Understanding the distribution and kinematics of ver- 
tically extended gas in disk galaxies is crucial for compre- 
hending the relevance for the evolution of disk -halo flows 
(e.g. the ga lactic fountain model presente d by Shapiro & 



|Field||1976| a nd refined byJB regman 1980; and the chim 
ney model of Norman & Ikeuchi|1989 1 , interactions with 



the intergalactic medium (1 GM) suc h as cold accretion 
of primordial gas ( Keres et al.||2005 Sancisi et al.||2008 1 
and interactions with neighboring galaxies. Furthermore, 
comprehension of the interactions between ISM energiza- 
tion from star for mation and accretion, i.e. feedback (e.g. 
Efstathiou||2000 ) is also imperative. 
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Galacti c gaseous halos have be en been observed in X- 
rays (e.g. TullmannetaL| 2006b[), radio continuum (e.g. 
Irwin et al.|1999[ ), d ust (e.g. Howk & Savage 1999]), ion- 
ized hydrogen (e.g. Rossa fc De ttmar 2003) as well as 
neutra l hydrogen (Hi) (e.g. S waters = et al.|l997 Poster 
loo et al.l|2007[) . For the radio continuum (e.g. ^ = 
et al.||200 6), diffuse ionize d gas (DIG) (e.g. |R< 
Dettmar 2003 ), X-ray (e.g. Tullmann et al i 
dust (Howk & Savage 1999D there is a correlation be- 
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,ossa Jz, 
2006al), and 



Despite the likely connection between several halo com- 
ponents and star formation, the prevalence and proper- 
ties of extra-planar H 1 in galaxies and any potential con- 
nection to star formation is not well understood. While 
there are many examples of H I shells and holes asso - 
ciated with star formation (e.g. Boomsma et al.l 2008), 
no correlation between the prevalence of widespread H 1 
halos and star formation activity has been established 
among multiple galaxies as has been done for other halo 
components. Clearly, this can only be remedied through 
greatly expanded observations. 

Another clue to halo origins lies in their kinematics. 
Trends concerning extra-planar H I, such as the presence 
and magnitude of any decreases in rotation speed with 
height (lags) or correlations with the kinematics of DIG 
layers are useful in determining its origins. One obser- 
vational issue is whether H 1 and DIG halos show the 
same kinematics, suggesting a common origin. Recent 
measurements hav e shown lags in multiple DIG layers 
leading to various models which at- 



(Heald et al. 



2007) 



tween the presence of these halo components and star 
formation in the disk, both locally, in regions through- 
out the galaxy as well as globally. This favors disk-halo 
flows akin to those in the aforementioned galactic foun- 
tain and chimney models as likely origins. 
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tempt to understand this gradient in terms of disk-halo 
flows and accretion of primordial gas. 
It has been shown that e ntirely ballisti c mode ls of disk- 



halo cycling, simulated in Collins et al. ( 2002 ) and Fra- 
ternali & Binneyl (I2006J) , are too simple in that they are 



unable to reproduce the observed kinematics. It is possi- 
ble the observed motions ma y partially be due to pressure 
gradients, magnetic tension ( Benjamin|2002 ) or external 
influences and f eedba ck such as those described in lFrater-1 
nali & Binney (2008). Additionally, "baroclinic" (where 
gas pressure does not depe nd on density alone) hydro- 



static models considered in Barnabe et al. (2006), have 



been able to reproduce the observed lag in MGC 891 
as has the hydrodynamic simulation of disk formation 
by Kaufmann et al. (2006). It is likely that halos form 
and evolve via a combination of disk-halo flows as well 
as external influences. 
Of further interest is the radial variation of lags. Basic 



Zschaechner et al. 



geometric arguments predict a shallowing of the l ag with 
increasing ra dius as may be seen in Figure 3 of Collins 



et al. (20021, but this could be complicated, for exam 
pie, by pressure gradients (Ben jam in 2002T). This will be 
discussed in greater detail in § |6.2| 

Detailed kinematic modeling of deep observations of 
H I emission has been done for only a small number of 
galaxies includ ing NGC 891 (Swaters et al.|1997l|Ooster- 
loo et al.|2007b, NGC 5746 (|Rand fc Benjamin|2008|) , and 
NGC 2403 ( |Fraternali et al. | |2002[ ). Through the West- 
erbork Hydrogen Accretiohin L Ucal GAlaxies (HALO- 
GAS) survey ( |Heald et aT"1|2011[ ), which targets 22 edge- 
on or moderately inclined spiral galaxies for deep 10 x 12 
hour observations using the Westerbork Synthesis Radio 
Telescope (WSRT), we aim to substantially increase this 
number in order to ascertain information concerning the 
origins of neutral extra-planar gas. Through deep obser- 
vations and kinematic modeling of these galaxies, we will 
establish whether there is any connection between extra- 
planar gas and star formation as well as the degree to 
which contributions from external origins are relevant. 

As part of the HALOGAS survey, we present observa- 
tions and models of NGC 4244, a nearby edge-on Scd 
galaxy, which exhibits a thickened H I layer upon in- 
spection of the zeroth-moment map (p resented here , Fig- 
ures [l][2]) as well as in previous work ( |011ing|1 996). The 
first goaf of the modeling is to determine the cause of the 
observed thickening, whether it is due to the presence of 
an H I halo, flare, warp along the line of sight, or any 
combination of these possibilities. Second, we constrain 
whether a lag exists. Finally, we will identify local fea- 
tures, such as possible en ergy injection sites, which will 



be discussed further in 5 6.4 This work will aid in es 



tablishing trends concerning the presence of halos and 
lags and any connection between disk-halo interactions 
and/or accretion. In the case of NGC 424 4, which has 
a low star formation rate (0.12 M yr _1 ; Heald et al. 
2011 1, we are especially concerned wi th the potential ex- 
lstence of a lag. If the trend found by Heald et al. (20071 
for DIG layers in galaxies with low star formation rate 
resulting in steeper lags holds true for H I, then a steep 
lag is expected. 

2. OBSERVATIONS AND DATA REDUCTION 

Here we provide a brief explanation of the observa- 
tions and data reduction process. A detailed description, 
which applies to all of the H ALOGAS galaxies i ncluding 
NGC 4244, may be found in |Heald et aT|p0TT ). 

Observations of NGC 4244 were primarily obtained as 
part of the ongoing HALOGAS survey. To best observe 
faint extended emission, the array was used in Maxishort 
configuration with baselines ranging from 36 m to 2.7 km 
for 6x12 hours. Of the 10 fixed antennas, 9 were used 
and were spaced at 144 m intervals on a regular grid. 
Addi tionally, 4 x 12 hours of archival data (Dahlem et al. 
2005 1 , observed in the traditional Westerbork variable 
configuration are used. The total bandwidth is 10 MHz 
with 1024 channels and two linear pol arizations. 

Data reduction was done in Miriad ( jSault et al. 1995). 
The H I data were imaged with a variety ot weighting 
scheme s, creating m ultiple data cubes. Clark deconvo- 
lution (Clark 1980) was performed using mask regions 
defined on the basis of unsmoothed versions of the data. 
Offline Hanning smoothing led to a final velocity resolu- 



tion of 4.12 km s . The la rms noise in a single channel 
of the full resolution cube is 0.22 mJy bin -1 , or corre- 
sponding to a column density of Njji = 3.5 x 10 18 cm -2 . 

For the primary modeling, we use a cube with a ro- 
bust parameter of for intermediate resolution and sen- 
sitivity, with a 21" x 13.5" beam, having a position an- 
gle of —0.7°. Additionally, to maximize sensitivity to 
faint extended emission we use a Gaussian uv taper cor- 
responding to 30" resolution in the image plane to cre- 
ate a second cube. Further smoothing to 60" resolution 
yielded negligible additional extended emission. Arcsec- 
onds are converted into parsecs via multiplying by 21.35 
pc arcsecond -1 [obtained by assuming a distance of 4.4 
Mpc (Tablell])] , resulting in a resolution of approximately 
450x290 pc for the full resolution cube. 

A summary of observational and data reduction pa- 
rameters may be found in Table [2] 

3. THE DATA 

The full resolution zeroth-moment map is produce d us 
ing the Gronigen Ima ge Processing System (GIPSY; 



van 



|der Hulst et al.|1992| task moments after applying a mask 
to the full resolution cube based on the 3<r limit of a 90" 



smoothed c ube. This is shown overlaid onto Ha (Hoopes 
et al |1999j ) and 24 ^m MIPS (NASA/IPAC Infrared Sci- 



ence Archive) images in Figure [T] as well as in grayscale 
in Figure [2J A warp component that is perpendicular 
to the line of sight is apparent. Also, as will be shown 
through modeling described in § |4j the minor axis extent 
is largely due to a warp component along the line of sight. 
Minor asymmetries may be noted, such as slight differ- 
ences in the warped disk in the approaching (northeast) 
and receding (southwest) halves as well as more radi- 
ally extended emission in the receding half. Additional 
localized asymmetries and irregularities are also present 
throughout. The 30" map (not shown) displays the same 
indications of a warp as well as the previously mentioned 
asymmetries. Even with the increased sensitivity of this 
map, it does not reveal any substantial faint emission not 
already seen in the full resolution cube. 

Full resolution representative channel maps are pro- 
vided in Figure [3j Note the faint streak-like emission 
in the lowest contours extending radially outward from 
the center of the galaxy in the 310 km s -1 and 322 km 
s -1 channels. Figure |4] shows a position- velocity diagram 
along the major axis (Tv diagram) with the rotation curve 
derived through modeling. The previously mentioned 
asymmetries are also easily seen in the channel maps as 
well as the lv diagram. 

Careful inspection of the data lays the foundation for 
the modeling process described below. 

4. MODELS 

To extract information from the data useful in con- 
straining the H I distribution and kinematics, including a 
possible lag, model data cubes are created using GIPSY. 
Using the task galmod, a model galaxy is divided into 
concentric rings where the column density, rotational ve- 
locity, scale height, velocity dispersion, inclination and 
position angle for each are specified. Additionally, we 
use a version of galmod, modified by one of us (G. H.) 
to allow a linear gradient in rotational velocity with 
height. These tasks assume axisymmetry, which is onl 
approximately true for NGC 4244 (best seen in Figures 
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Fig. 1. 
by 43°. 
The Ha 



The Hi zeroth-moment map overlaid on Ha (top) and MIPS 24 fim (bottom) images. Each image is rotated counter-clockwise 



Hi contours begin at 6.4 X 10 



image is supplied by Rene Walterbos from|Hoopes et al.|(1999) 



and increase by factors of 2. The Hi beam is shown in the lower left-hand corner in white. 



TABLE 1 
NGC 4244 Parameters 



Parameter 



Value 



Reference 



Distance (Mpc) 


4.4 a 


Systemic velocity (km s~ 1 ) 
Inclination 


211 
88° 


SFR (M yr" 1 ) 
Morphological Type 
Kinematic Center a (J2000.0) 
Kinematic Center <5 (J2000.0) 


0.12 

Scd 

12hl7m29.90s 

37d48m29.00s 


Optical Radius (kpc) 

Total Atomic Gas Mass (1O 9 M ) 


10.4 b 
2.5 


Total Atomic Gas Mass (10 9 M Q ) 


2.2 



H eald e t al. (2011) 



Oiling (1996), This work 
This work 
Heald et al. (2011 



|de Vaucouleurs et al.| ( 11991 
i ms work 
This work 



Fry et al. 
This 



1 1999 I 



Oiling (1996) 



a Distance is the median value of distances found on the NED database, excluding those 
obtained using the Tully-Fisher relation. 

This value is adjusted for the difference in distance of 4.4 Mpc used in this paper and 
3.6 Mpc used in |Fry et al~1jl999[ 
c Includes neutral He via a multiplying factor of 1.36. 

Scaled to our assumed distance. 



through El). Therefore, the approaching and receding 
halves are modeled separately, which greatly improves 
the models for each half but does not completely elimi- 
nate issues due to asymmetries as will be demonstrated. 
These remaining issues however, do not impact our final 
result concerning the existence and magnitude of a lag. 
For initial estimates of the column densities as well as 



the rotation curve, we used the models of Oiling (1996), 
who found evidence for a flaring gas layer. Given the mi- 
proved depth of our data by comparison (with the rms 
noise for this data being 0.22 mJy bm _1 compared to 1.9 
mJy bm _1 for a cube of similar resolution (28.3x10.1") 
used by Oiling) as well as differences in modeling tech- 



niques used, it is possible for the results to differ. 

Rough initial estimates for the kinematic center and 
systemic velocity were found using the zeroth and first 
moment maps. Estimates for the position angles were de- 
termined via inspection of the zeroth-moment map and 
were immediately refined prior to altering additional pa- 
rameters. The estimated column densities are refined 
first by matching both the shape and total flux in the 
vertical profile (Figure pi made by including emission 
only within a major-axis distance of 5' from the center 
to avoid effects such as the warp and tapering of the disk 
found at larger radii) and zeroth-moment maps of the 
model with those of the data. 
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18 m 17 m 20 s 12 n 16 m 40 B 

H.A- (3000.0) 

Fig. 2. — A zeroth-moment map of NGC 4244 showing the slice 
locations for the bv plots in Figure m\ (solid) , as well as the range 
of the evenly spaced slices in FigurcTnol (d ash ed). Bo xes a and b 
correspond to the regions shown in Figures |ll| and |13| respectively. 
Contours begin at 6.4 X 10 19 cm -2 and increase by factors of 2. 
The H I beam is shown in the lower left-hand corner in white. An 
encircled dot and cross denote the approaching and receding halves 
respectively. 



TABLE 2 
Observational and Instrumental Parameters 



Parameter 



Value 



Observation Dates - HALOGAS 



Observation Dates — Dahlem 



Pointing Center 

Number of channels 
Velocity Resolution 
Full Resolution Beam Size 

RMS Noise - 1 Channel (21x13.5") 
RMS Noise - 1 Channel (30x30") 



2009 Jul 13 

2009 Jul 16 

2009 Jul 18/Aug 1 

2009 Nov 15 

2009 Nov 27 

2003 Feb 13-14 

2003 Feb 19-20 

2003 Feb 20-21 

2003 Apr 13-14 

12hl7m29.901s 

37d48m29.00s 

68 

4.12 km s _1 

21x13.5" 

450x290 pc 

0.22 mjy bm" 1 

0.37 mjy bm _1 



Provided a good fit for previous steps, the rotation 
curve is then honed primarily through the use of lv di- 
agrams, tracing emission along the edge with the high- 
est velocities (terminal side). These, along with minor 
axis (bv) position- velocity diagrams, as well as channel 
maps of models and data must all be examined to fine 
tune the above parameters in the models, while at the 
same time ensuring that the vertical profiles and zeroth- 
moment maps are well matched. Although all plots are 
considered throughout the process, the bv diagrams and 
channel maps are the most heavily relied upon after the 
initial modeling stages. 

As will be demonstrated, the key to this method is 



to start with the simplest possible model and then add 
features such as warps and flares - one at a time - to 
fully understand the individual contribution of each ad- 
ditional degree of complexity. This allows for subtle dis- 
tinctions between features, which otherwise may be mis- 
interpreted. In the end, an acceptable model is one that 
fits all of the plots well. The quoted errors indicate visu- 
ally estimated uncertainties rather than formal la error 
bars. 

It should be noted that ne w, semi-automate d tilted 
ring fitting software (TiRiFiC) ( Jozsa et al.|2007 ), which 
allows for a chi-square fit of models to observed data 
cubes is now available. Advanced models for NGC 4244 
were examined in both TiRiFiC and galmod and each 
produced comparable results. Given the advanced stage 
of modeling already completed via galmod, as well as 
the localized features in NGC 4244, TiRiFiC was not 
used extensively in this case. However, the advent of 
TiRiFiC's automated processes will help to expedite the 
fitting of the remaining galaxies in the HALOGAS sam- 
ple. 

4.1. Features Considered While Modeling 

The bv diagrams shown in Figure [6] and the channel 
maps in Figure [7] illustrate a number of the features that 
were noticeably sensitive to the model parameters we en- 
deavor to match. There is a notable "T" shape in the 
panels closest to the galactic center in the bv plots in 
Figure [6j which morphs into a "V" shape in the outer 
panels (schematically indicated in certain panels of Fig- 
ure [6]) . Additionally, the contours on the systemic side 
are Flattened. There is also a substantial lopsidedness 
best seen at —5.1' and —6.9' in Figure [6] around the mid- 
plane, primarily due to the component of the warp per- 
pendicular to the line of sight, but also partially due to 
asymmetries. 

The channel maps displayed in Figure [jj reiterate the 
need for a warp perpendicular to the line ofsight. Promi- 
nent features that were sensitive to model parameters 
considered during the modeling process include the spac- 
ing and angle of the contours at the tips of each plot 
closest to the center of the galaxy. Also considered is the 
approximately 45 degree slant along the outermost edges 
at more extreme velocities and at mostly negative minor 
axis offsets in the approaching half and positive offsets 
in the receding half. 

In both Figures [6] and [7J the visual assessment of the 
quality of the fit between the models and data is indi- 
cated in each panel via symbols explained in the caption 
of Figure [6j The assignment of each symbol heavily relies 
on the criteria described above. 

Some quantities are kept constant throughout all of 
the models for the entire process. These are the systemic 
velocity (244 km s _1 ), kinematic center (12hl7m29.90s, 
37d48m29.00s) and the run of the position angle derived 
from the component of the warp perpendicular to the 
line of sight, column density profile, rotation curve (both 
shown in Figure p| and velocity dispersion (12 km s _1 
decreasing to 10 km s _1 at a radius of 7.7' or 10 kpc). 

4.2. Individual Models 

Several possible models are considered. Some are elim- 
inated quickly upon inspection of the vertical profile or 
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Fig. 3. — Representative channel maps of the data with a resolution of 21" X 13.5" . Contours begin at 3cr (0.66 mjy bm x ) and increase by 
factors of 3. Negative contours are also shown with dashed lines. The velocities in km s — 1 are given in each panel and Sv is approximately 
3 resolution elements. Note the slanted outer edges present in both halves, which correspond to the projection of the warp perpendicular 
to the line of sight. Arrows indicate the pronged streak discussed in the text in the lowest contours of the channels corresponding to 310 
and 322 km s — . 
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Fig. 4. — Major axis position-velocity diagram. Contours begin 
at 2o~ (0.44 mjy bin -1 ) and increase by factors of 2. Asymmetries 
are once again prevalent in the two halves. The rotation curve for 
each half is shown with crosses. Arrows near the horizontal axis 
indicate the approxi mate st arting point of the radial decrease in 
the lag discussed in § |4.2.6| The warp component across the line 
of sight creates the appearance that the rotation curve extends 
beyond the data, but at heights above the midplane, these points 
are constrained. 
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Oiling (1996) 




-3 -3 -1 +1 +3 +3 

Offset R.A. (ARCMIN) 

Fig. 5. — Log-scale vertical profiles of the data, flaring model and 
Olling's 1996 model with a constant inclination for both halves. 
The other models in this paper are not shown as they match the 
data closely, to within one line width. Notice how the data are fit 
with a single exponential and do not display prominent emission 
in the wings. 



zcroth-momcnt map, while others require substantially 
more effort to discern between them. Vertical distribu- 
tions are exponentials in all models. Exponential scale 
heights vary among models and are listed in Table [3| All 
of the models, unless otherwise noted, match the zeroth- 
moment map and vertical profile reasonably well, as will 
be graphically shown in Figure [9] and in the following 
discussion. 

4.2.1. One-component Model 

The simplest model involves a single disk with a slight 
warp perpendicular to the line of sight [deviating by 4- 
5 ° from the inner disk (Figure [8]) and beginning at a 
radius of 12 kpc]. For such a model, all other quantities 
aside from column densities and rotational velocities are 



TABLE 3 

Model exponential scale heights 

rounded to the nearest 25 pc. 



Parameter 



Value 



scale height 
scale height 
scale height 
scale height 
scale height 
scale height 



1C (pc) 
F(pc) 
W (pc) 
WF (pc) 
2C (pc) 
2CW (pc) 



600 

550 -> 1100 a 

550, 575 b 

550 -> 1100 a 

525, 1050 c 

425, 850 c 



a the first value applies to radii within the 
optical radius. Beyond this radius, the scale 
height linearly increases until the final value 
(twice the initial) is reached. 

Values for the approaching and receding 
halves respectively. 

c Values for the thin and thick components 
respectively. These are consistent for both 
halves. 

held constant with radius. The inclination is found to 
be 88 ° for the best fit to the vertical profile. For the fi- 
nal model, the scale height is about 600 pc. This simple 
model (referred to as "1C" in figures) cannot be imme- 
diately eliminated based on the fit to the vertical profile 
and zeroth-moment map (Figure p)l) . However, one may 
note that in the bv plots in Figurejol the one-component 
model lacks sufficient flux on the systemic relative to the 
terminal side in each panel. As can be seen in the chan- 
nel maps in Figure [71 shown most prominently in the 
plots closest to the systemic velocity and noted with an 
arrow, the edge of the diagram away from the center of 
the galaxy comes to a definitive point at a positive minor 
axis offset on the approaching half, and a negative offset 
on the receding, which is not seen in the data. No other 
combination of inclination and scale height improves this. 

4.2.2. Warp Along the Line of Sight 

Adding a projection of the warp along the line of sight 
(referred to as "W" in figures), by gradually decreas- 
ing the inclination in the outermost rings (shown in Fig- 
ure 8l), matches the observed "T" shape from the data 
in the central bv diagrams in Figure [6] significantly bet- 
ter than the 1C model. This is the result of bringing 
flux from outer rings with lower projected velocity into 
higher latitude regions. If the inclination decreases too 
quickly, the "T" shape becomes excessive. 

In contrast to the one-component model, the channel 
maps in Figure [7] show the flux is pushed upward and 
at high z pushedaway from the center in the major axis 
direction, in panels with velocities close to the systemic 
velocity in the warp model, which matches the data more 
closely. However, if too much of a warp is present in 
the model, then analogous to the "T" shape verging on 
becoming too "Y"-shaped in the bv diagrams, there will 
be an indentation on this edge. This is not seen in the 
data. 

The best fit LOS warp component deviates by 8 degrees 
from the inner disk and as in the PA warp, starts at a 
radius of 12 kpc (Figure [8]). 

4.2.3. Flaring Model 

Upon initial inspection of the zeroth-moment maps, it 
may be noted that the edges in the outer radii taper 
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Fig. 6. — Minor axis position- velocity diagrams of NGC 4244 and models. Slice locations are those displayed by solid lines in Figure [2] 
and velocity ranges vary from panel to panel. Contours begin at 3cr (0.66 mjy bm _1 ) and increase by factors of 3. All models include a 
projection of a warp across the line of sight. The W model is favored and others are shown for illustrative purposes only. The symbols 
provided in the lower right-hand corners of selected panels indicate quality of fit to the data judged by eye. Crosses indicate a very poor 
fit, blank a somewhat poor fit, boxes a reasonable one, and asterisks a good fit. The "V" shape and flattened systemic side are emphasized 
with gray dotted lines, while the "T" shape is shown in black and white lines. (A) displays models with no lags, while (B) shows analogous 
plots, but with the addition of lags as well as the omission of the 2C model. The magnitudes of lags are —7 km s — 1 kpc - 1 for both 
halves of 1CL as well as FL, —9 km s _1 kpc -1 for the WL and WFL models, and —9 km s — 1 kpc -1 for the 2CWL model. Due to local 
asymmetries, the emission extending on the negative minor axis offset near 125 km s — * in the panel corresponding to 6.9' is not considered 
while modeling. Additionally, on the receding half, the side with a positive offset is given less weight in the fitting process for the same 
reason. 



significantly in minor axis extent. Because of this, a sub- 
stantial flare can be ruled out immediately as it would 
cause the zeroth-moment map to become too boxy and 
fail to fit this aspect of the data. However, this is not to 
say for certain that a modest flare may not be present. 
For this reason, a representative flaring model (referred 
to as "F" in figures) is created by starting with the one- 
component model and increasing the scale height out- 
ward radially. This increase begins at the optical radius 
(10.4 kpc; Tablell]) and increases linearly, approximately 
doubling by a radius of 15 kpc. Comparison of zeroth- 
moment maps from the data (Figure Kn and model in- 
dicates that a modest flare yields a slight improvement 
compared to the one-component model. Unlike the warp 
model, the modest flare does not eliminate the rounding 
on the systemic side of the bv plots (Figure]!]) or the thin- 
ning on the outer edges of the channel maps seen in the 
one-component model. A more extreme and less physi- 



cal flare (abruptly doubling the scale height at the optical 
radius and then holding it constant at larger radii) helps 
to reduce these problems in the bv plots and channel 
maps, but once again, causes the zeroth-moment map to 
become boxy and an obvious mismatch to the data (not 
shown) . 

The pure flare model clearly does not match the data 
in the bv plots. However, adding the same modest flare 
to the warp model (referred to as "WF" in figures) does 
not appear to have a substantial effect; in fact, it is al- 
most identical to the warp model as shown in the zeroth- 
moment map in Figure [9] as well as the bv plots in Fig- 
ure [6J This indicates that such a flare could be present in 
NGU4244, but that the component of the warp along the 
line of sight is dominant and effectively hides most of its 
contribution. Since we are unable to determine whether 
or not the flare exists due to the minimal contribution 
it would have compared to the warp, and for the sake of 
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clarity, it is omitted from further models. 

4.2.4. Two-component Model 

H I halos have been observed in several galaxies, with 
NGC 89f among the most notable (Oosterloo et al. 2007). 
One of the primary goals of HALGGAS is to determine 
the frequency of occurrence and origins of such halos. We 
therefore examine the possibility of such a halo in NGC 
4244. 

The addition of a halo to the model is made via com- 
bining a second, thicker exponential component to the 
first, keeping the inclination of both at 88° (referred to 
as "2C" in figures). Upon inspection of the vertical pro- 
file and zeroth-moment map, there is no immediate mo- 
tivation for a second component, as reasonable fits can 
be made with less complicated models. For both halves 
of NGC 4244, multiple ratios of scale height and flux in 
each component were tested in attempts to match the 
vertical profile, but none was successful. In fact, when 
considering only the vertical profile, a two-component 
model cannot be made to match the data, as during the 
modeling process, the components tend to converge to a 
single scale height, resulting in a one-component model. 

However, since the existence of a halo is of primary 
concern, we have made substantial attempts to optimize 
this model. Since the vertical profile and zeroth-moment 
map provide no useful constraints for the relative scale 
heights of the two-components or for the percentage of 



flux contained within the halo, we use illustrative fixed 
values, with the scale height of the extended component 
twice that of the thin component, and with 25 percent 
of the flux residing in the extended component. These 
values are chosen because they represent a reasonable 
scenario compa red to other galaxies with halos ( |Ooster- 
looet al.|[2007| . 

Through optimization of this two-component model, it 
is found that a reasonable fit to the data is impossible 
without including a warp component along the line of 
sight. A halo will only add to the complexity of the model 
and if a warped model without a halo is sufficient (or in 
this case, slightly better) , then the simpler model should 
be used. Thus, we reject the two-component model with 
or without such a warp. We do however, retain a two- 
component model with a warp component along the line 
of sight ("2CW") in figures throughout this paper for 
illustrative purposes only. 

4.2.5. Improvement With the Addition of a Lag 

At this point, prior to the addition of a lag, we place 
our confidence in the single component LOS warp model. 
However, as may be seen in Figures [6] and [7[ the opti- 
mized models with no lag do not exactly match the data. 
Adding a lag at this stage will best illustrate its contri- 
bution to the models. We continue to include the one- 
component model, the flaring model with a warp compo- 
nent along the line of sight, as well as the two-component 
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Fig. 7. — Channel maps of NGC 4244 compared with models. Contours begin at 3cr (0.66 mjy bm~ 1 ) and increase by factors of 3. 
Velocities are given in the first column in kilometers per second. All models include a projection of a warp perpendicular to the line of 
sight. Lag models are denoted with an L. (See Figure J6] caption for the lag magnitudes.) The pure flare and pure two-component models 
are not included as they are shown to be poor matches in the bv plots. Finally, the WF model has also been excluded as it is similar to 
the warp only model. The definition of symbols follows that set forth in the caption of Figure [m The approaching half is shown in (A) 
while the receding half is shown in (B). 



model with a warp component along the line of sight for 
illustrative purposes as well as to demonstrate the relia- 
bility of the result concerning the lag. 

The columns in Figure [6] denoted with "L's" in the ti- 
tles show the models as before, but with optimal lags for 
each model. Upon inspection of these plots, it becomes 
evident that the addition of a lag is an almost univer- 
sal improvement. By adding a lag, flux in the models is 
drawn to the systemic side, thus duplicating the afore- 
mentioned "T" and "V" shapes in the data. Further- 
more, the edge of the systemic side becomes flatter in 
almost all of the panels as in the data, whereas in the 
no lag models, only the models including a line of sight 
warp component appear relatively flat. Additionally, the 
channel maps are improved. They narrow in the z di- 
rection at velocities further from systemic, more closely 
matching the data. 

The warped model with the addition of a lag of 
km s _1 kpc^ 1 and — 9±2 km s _1 kpc -1 in the approach- 
ing and receding haves respectively is a reasonable fit. 
Error bars are estimated via visual inspection of the 



q+3 

y -2 



range of lags which could potentially match the data and 
under the (inexact) assumption that the models have no 
additional uncertainties (e. g. uncertainty in column 
density or rotational velocity), as making such estimates 
would be unwieldy. 

For the two-component warp model, a lagging halo is 
simulated by adding a lag in only the thick component. 
The addition of the lag substantially improves the model 
for the receding half with a lag as steep as —9 km s _1 
kpc -1 , but fails to adequately duplicate the "V" shape 
in Figure [6] in the approaching half. Even with a lag 
of as steep as —28 km s _1 kpc -1 this shape cannot be 
duplicated. Given this difficulty, a lag of —9 km s _1 
kpc -1 is used in the final model to optimize the fit for 
the overall shape of the data. 

The improvement due to the lag is present in all mod- 
els. This attests to the robustness of the lag in that, 
regardless of the optimal morphological model, there are 
features which cannot be re-created without it. Once 
again, the one-component warp model with the addition 
of a lag appears to be the best match to the data. 
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Fig. 8. — Shown are parameters used for the models. The ro- 
tation curves, column densities and position angles presented are 
consistent for all models. The inclinations are those used for the 
warp and two-component models. 



4.2.6. The Radial Variation of the Lag 

Upon establishing the existence of a lag, it must be 
characterized by not only the magnitude of the gradient 
in the z direction but also its radial variation. Using the 



one-component warped model, bv diagrams further from 
the center of the galaxy th an those shown in Figure 
are examined in Figure [10| using a 30x30" cube, wit 



the range of slice locations corresponding to the dashed 
lines in Figure [2| The data, at these large radii, are 
best fit with a shallowing lag, decreasing in magnitude 
to — 5±2 km s _1 kpc" 1 and — 4±2 km s _1 kpc -1 near 
a radius of 8' (10 kpc) in the approaching and receding 
halves respectively. This is most easily seen on the ter- 
minal side of each panel: the data contours become more 
rounded, whereas the model with a constant lag retains 
too much of a "V" shape, with data contours extend- 
ing further on the terminal side than those of the model. 
However, on the negative offset side of the receding half, 
the shallowing of the lag may be even more substantial. 
This is best seen in panels corresponding to -8.2' and -9'. 
Uncertainties are still quantified based on the positive 
offset side, and —4 km s" 1 kpc" 1 is only given as an up- 
per limit for the side with negative offset. There is no 
strong indication of any radial variation at radii smaller 
than approximately 6' (7.5 kpc). It should be noted that 
the S/N of the data is considerably diminished at and 
beyond a radius of 10' and determining the lag this far 
from the center approaches the limits of our modeling 
capabilities. The error estimates are once again by eye. 

5. DISTINCT FEATURES 
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Fig. 9. — Zeroth-moment maps of the data and models. The 
model plots are comprised of both the approaching (positive offset) 
and receding (negative) separately modeled halves. Contours begin 
at 6.4 X 10 19 cm -2 and increase by factors of 2. 

Aside from the global characterization of H I kinemat- 
ics and distribution, it is a goal of HALOGAS to charac- 
terize localized features which may be accretion events or 
indications of disk-halo cycling. Furthermore, there exist 
prominent features in the data we wish to demonstrate 
are not global in the sense that our models do not at- 
tempt to reproduce them. Perhaps the most significant is 
a shell-like feature located in the western quadrant of the 
approaching half. An arcing feature is seen in a zeroth- 
moment map made by including channels ranging from 
149 km s _1 to 175 km s _1 , extendi ng f rom roughly 3.6' 
to 5.1' along the major axis (Figure 111. Corresponding 
extended emissi on is also observed as line splittin g in t he 



6.4 



bv plots (Figure [12| , which will be discussed in 

An additional note wor thy feature is observed in the 
receding half (Figure 13), roughly 4.3' to 6.8' radially 



from the center, and present in channels corresponding to 
281-302 km s^ 1 . This feature displays an elongated path 
curving away from the major axis in the Hi contours, 
starting from the midplane and extending beyond 1.5' 
above the midplane. Furthermore, it appears that there 
may be a hole or vent in the H I at this location, which 
may correlate with Ha features (Fig ure [T]). Details of 



6.4 



this feature will be discussed in 

A narrow, pronged, faint streak may be seen from 302 
km s _1 to 355 km s" 1 in the channel maps of the receding 
half (Figure J3J). Beyond 355 km s~ x it becomes difficult 
to discern the streak from the regular emission of the 
galaxy which can be modeled. However, by comparing 
the H I along the mid-plane in these channels for each 



half, it is seen that the receding half extends only slightly 
further than the approaching. Thus, it is possible that, 
rather than there existing a streak of additional H I, there 
may instead be a dearth of H I at moderately high z in 
these channels. This feature will not be discussed further 
as we have no clear explanation. 

Several small-scale shell-like features are also seen in 
both halves of NGC 4244, but are not shown or discussed 
further as they do not have substantial effects on the 
modeling and appear to have no immediate connection 
to features in the plane of the disk. 

6. DISCUSSION 

6.1. The Warp in NGC 4244 

Our models show that a warp starting around a ra- 
dius of 12 kpc exists in NGC 4244. This warp is not only 
seen in the zeroth-moment map, forming a classic shallow 
integral shape, but also along the line of sight through 
our more in-depth modeling. These two-components of 
a warp, while representing tilts around axes off-set by 
90° from each other, are similar in amplitude (8° and 
4-5 ° ) and starting point (Figure [8]) . Considering both 
contributions, the angle included Hoy the spin vector of 
the inner, flat disk and the spin vector at the outermost 
radius, 17.5 kpc, for which we can reliably determine in- 
clination and position angle (i. e. the warp amplitude) 
is 9-9.5°. The starting point, being near the optical ra- 
dius is consistent with t he first two r ules of commonly 



observed warp behavior (Briggs 1990). This agreement 
between two independently modeled quantities as well 
as the consistency with past observations attests to the 
reliability of the models. 

6.2. Impact on Lag Trends 

In addition to the presence of extra-planar gas, trends 
concerning its properties such as the existence and mag- 
nitude of lags are of use in determining its origins. Al- 
though we do not detect a vertically extended H I halo 
in NGC 4244, we do detect a lag in a thickened H I layer 
having a scale height of approximately 550-575 pc. This 
indicates that substantial extra-planar gas need not be 
present in order for a lag to be observed, as is potentially 
seen in UGC 1281 ( |Kamphuis et al.||2011[ ), although for 
that galaxy, a more substantial line of sight warp com- 
ponent is favored in place of a lag. 

Any correlations with DIG la yers and their kinemat- 
ics are also critical. Results from Fraternali et al. ( 2004 1 
show agreement between the kinematics of the Hi and 
DIG layers off the plane in NGC 2403, provid ing an 
early indica tion that they may be connected. |Heald| 
et al. ( 2007 ) discovered a trend [albeit for only 3 galax- 
ies: NGC 891 (-17.5 km s" 1 kpc" 1 ), NGC 4302 (-30 
km s- 1 kpc" 1 ) and NGC 5775 (-8 km s" 1 kpc" 1 )] of a 
shallower lag with increasing star formation activity for 
extraplanar DIG. At this time, Ha data for NGC 4244 
is being analyzed by additional HALOGAS members to 
determine the DIG kinematics and results are forthcom- 
ing (Wu et al. 2011, in prep). Although we cannot say 
for sure, the low star formation rate of 0.12 M yr _1 
in NGC 4244, which is substantially lower than those in 
the galaxies previously listed, indicates that if a lag in 
the DIG exists, then it should be steep compared to the 
trending galaxies. However, the lag in H I in NGC 4244 
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Fig. 10. — bv diagrams showing the warp model with no lag (left-most column), a best estimate for a constant lag (center), and a best 
estimate for a radially varying lag (right). The cube represented here has a resolution of 30". Contours begin at 3<r (1.1 mjy bm" 1 ) and 
increase by factors of 3. Data contours are black and model contours are superimposed in white. The lags in km s — 1 kpc - 1 for each panel 
of the varied warp model are given, while the lag for the constant warp model is given in the top panel only. The offset of the slices from 
the center in arcminutes is given in the left-most panels and correspond to the dashed lines in Figure [2] The quality of fit of individual 
panels to the data is indicated by the convention initiated in Figurem] Boxes are placed to show the more subtle differences between panels. 
Differences which are more readily apparent to the eye are not explicitly marked in this way. The primary issues to note are the slope of 
the model contours compared with those of the data on the terminal side as well as any shift up or down between them. Also considered 
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Fig. 11. — A zeroth moment map of the shell-like feature on the 
approaching half centered at 12hl7m42s, 37d53m9s and described 
in the text (top) and plotted over the Ha image (bottom) with the 
H I beam displayed in the upper right-hand corner of the top plot. 
The H I contours start at 4.7 X 10 18 cm — 2 and increase by factors of 
V2. The velocity range included is 162-174 km s _1 The box in the 
lower panel outlines the star forming region described throughout 
the text. 




Minor Axis Offset (ARCSEC) 

Fig. 12. — bv plot s at the location of the shell-like feature dis- 
played in Figure |11| The dashed lines in the center panel show 
line splitting at the center of the presumed shell, while the left and 
right panels correspond to the approximate locations of the shell 
walls. Contours begin at 0.56 mjy bm _1 and increase by factors 
of \/L8. 



is shallower than the DIG lags for NGC 4302 and NGC 
891 and slightly steeper than that found in NGC 5775. 
Given this, either the trend of low star formation with 
steeper DIG lags will be broken if the H I and DIG lags 
match, or the magnitudes will differ, possibly indicating 
d ifferent origin s. 



Heald et al. (2007) also found that, although the lags 
differed substantially in the three DIG layers studied, 
they are in better agreement with each other when ex- 
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Fig. 13. — A zeroth moment map of the curved feature on the 
receding half described in the text (top) , plotted over the Ha image 
(middle) as well as a zoomed-in image (bottom) of the boxed area 
in the Ha image with an enhanced grayscale showing potential 
outflow from the Hn region (in the boxed areas), which falls in 
a gap in the Hi. The Hi contours start at 3.7 X 10 19 cm~ 2 and 
increase by factors of \/2. The velocity range included is 281-302 
km s _1 . 
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Fig. 1. — An Hi channel map at 298 km s _1 . 3 and 7a Ha 
contours are overlaid on an H I grayscale to show the outline of the 
faint Ha emission. Notice the gap in the H I near the location of 
the potential Ha outflow also shown in Figure |l3] 
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pressed in terms of km s _1 per DIG scale height with 
an average value of —20 km s -1 per scale height. If one 
considers the NGC 4244 lags in such terms, the Hi lag 
found in NGC 4244 is dissimilar to those found in the 
DIG layers of previously modeled galaxies. Assuming 
the warp model, with scale heights of 550 pc and 575 pc 
in the approaching and receding halves respectively, the 
lag is roughly —5 km s~ x per scale height. 

However, it is useful to note that possibly differing 
magnitudes in the H I and DIG lags do not necessar- 
ily mean that they do in fact have different origins. The 
H I and DIG have substantially different distributions, 
with the H I having significantly greater radial extent as 
may be seen in Figure [l] Since the magnitude of a lag 
from an internal source is somewhat dependent on where 
the gas is initially launched, at least in ballist ic models 
HCollins et al.j|2002[ [Fraternali fe Binneyl2006) , it is safe 
to say that different radial density profiles, as are com- 
monly seen in H I and DIG layers, can potentially result 
in differences in the lags. 

We now compare the magnitude of the lags in NGC 
4244 with lags found in the H I halos of other galaxies. 
Gradients in Hi laye rs have been modeled for NGC 891 



[-15 km s- 1 kpc- 1 dOosterloo et al. 2007)LUGC 7321 



-25 km s" 1 kp^ 1 dMatthews fc Wood]2 003)l and the 
Milky Way [-22 km s l kpc" 1 dLevine et al|2 008), -15 



km s" 1 kp^ 1 (Marasco & Fraternali 2010&J. The lag 



measured for NGC 4244 is substantially shallower than 
these. The star formation rate in NGC 4244 of 0.12 M 
yr _1 is su bstantially lower than that of NGC 891 [3.8 
M Q yr" 1 dOosterloo et al.||2007|)1 and the Milky Way 
[0.68-1.45 M yr- 1 dRobitaille fc Whitney||2010|l, but 
comp arable to UGC 7321 [0.15 M yr" 1 QHeald et al 
2011 1]. Given the small number of galaxies mentioned 
above and the lack of any pattern concerning star forma- 
tion rate and the magnitude of a lag in extraplanar H I, 
it is still too early to establish a significant trend. 

6.2.1. The Relation Between the Lag and Star Formation 
Within NGC 4244 

Although the lag in the H 1 layer of NGC 424 4 does 
not fit the t rend with star formation rate seen by |Heald 
et al. (2007) for DIG layers, we can ask whether such a 
trend is apparent within NGC 4244: the lag is compara- 
ble on the two sides of the gala xy, but is this also true of 
the star formation rate? If the Heald et al.| ( |2007[ ) trend 
continues, then a large contrast in star formation rate 
would be unexpected. We consider the star formation in 
each half separately by using Ha (Hoopes et al. 1999) 
and 24 ^m MIPS data (NASA/IPAC Infrared Science 
A rchive). Using these, w e apply the method described 



Kcnnicutt et al. 



( ]2009| ) to correct for extinction. We 
estimate the uncorrected Ha flux in each half by mea- 
suring the percentage of the total Ha flux in each and 
th en using the total Ha luminosity for NGC 4244 found 
'Kennicutt et al.| (|2008[) of 10 40 erg s _1 . Accounting 



m 
for 



the slight difference m distances used, these estimates 
are 5.4 x 10 39 and 4.4 x 10 39 erg s" 1 in the approach- 
ing and receding halves respectively. (This ordering will 
remain consistent throughout this section.) Assuming a 
monochromatic flux (L = XL\) for the 24, 70 and 160 
/zm data, the total infrared (TIR) fluxes are found to 
be 1.1 x 10 42 and 1.0 x 10 42 erg s" 1 . By assuming a 



Salpcter IMF and then applying a conversion between 
line fluxes and st ar formation rates supplied by |Kenni-| 
cutt et al. (20071, the star formation rate for each half 
may be determined: 

(1) 



SFR 



M 



©yr" 



9x 10-^[L na + aL Tm ] 



(where a = 0.0024 and units are in erg s 1 ). These are 
found to be 0.064 and 0.054 M Q yr -1 . The star forma- 
tion rate in the approaching half is approximately 15% 
higher than that in the receding half. However, the un- 
certainties in both the lag magnitude and the local star 
formation rates do not allow for an unambiguous co nclu- 
sion regarding the trends seen by Heald et al. (2007), but 
there is no indication of deviation from it. 

6.2.2. The Radial Variation of the Lag 

Now we consider the radially outward decrease in the 
magnitude of the lag, which has yet to be fully explained. 
This decrease begins just within the optical radius (near 
7' or 9 kpc) and continues until just past the start of the 
warp. Thus, it may be concluded that this is an effect 
independent of the warp. By comparing the span of this 
radial variation in lag with the distribution of star for- 
mation as traced in Ha and 24 /jm emission (Figure [l]), 
The onset of the decreasing lag magnitude coincides with 
a sharp decline in star formation, again going against 
th e tentative tr e nd fo r extraplanar DIG rotation found 
((20071. 



by Heald et al. 
sume that star 



However, it is premature to as- 



ormation is the only process affecting the 
magnitude of the lag and other factors may be the cause. 
Simple geometric arguments predict a shallowing of the 
lag at large radii. The most basic interpretation of the 
geometry, without any consideration of gas dynamical 
effects, considers the relationship between the rotational 
velocity and the gravitational potential: gas rotating at 
high z is further from the gravitational center than the 
gas at the same radius in the midplane, resulting in a 
slower rotational velocity. This effect is most noticeable 
at smaller radii, where a certain height above the disk is 
a larger fraction of the radial distance from the center, 
resulting in a larger relative change in the gravitational 
potential, and thus a steeper lag. 

Additionally, a radial pressure gradient c ould affect 
the la g and its radial variation. According to Benjamin 
( 2002 1 , a radial pressure gradient directed inward may 
cause a steepening of the lag. In contrast, an outward 
pressure gradient may result in a shallowing of the lag 
or even an increase in rotation velocity. Such gradients 
could also vary with radius. Given the uncertainty in the 
pressure gradients in question, as well as the potentially 
substantial impacts even small pressure gradients may 
have, it is difficult to say to what degree these will alter 
the magnitude of any given lag. 

Observa tionally, a radially var ying lag is also seen in 

NGC 891 dOosterloo et al.||2007|), where the gradient in 

r 



the inner regions is 
-14 km s -1 



-43 km s 1 kpc and decreases to 



km s 1 kpc" 



kpc 1 in outer radii. This shallowing of 2.5 
2 is comparable to what the models indicate 



for NGC 4244, although the uncertainties (as judged by 
eye) make this a difficult comparison. Furthermore, the 
shallowing of the l ag in NGC 891 was se en closer to the 
center in both Ha dKamphuis et al.|2007 ) and H 1. Thus, 
the origins of the shallowing of the lags in NGC 4244 and 
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NGC 891 may be different. 

Also displaying radial variati on in the lag is the DIG 
in one quadrant of NGC 4302 ( |Heald et aL]|2007[ ), but 
the trend found in that galaxy was opposite from what 
we detect in the H i in NGC 4244. NGC 4302 displayed 
a substantial steepening of the lag, going from approxi- 
mately —23 km s _1 kpc -1 to almost —60 km s _1 kpc -1 
at radii outward of 4.25 kpc. 

Ultimately, a larger sample of galaxies where this phe- 
nomenon is observed and modeled will greatly aid in de- 
termining its cause and nature. 

6.3. Impact on H I Halo Trends 

There is substantial evidence against the existence of 
an extended H I halo in NGC 4244, most notably the 
failure of the addition of a second vertical component by 
itself or accompanied by a warp along the line of sight 
to improve the fit of the vertical profile (Figure [5| , and 
that the appearance of Figure [l] does not differ substan- 
tially from the 10 x shallower Oiling data. Furthermore, 
a second component added to the warp model yields no 
overall improvement to the bv plots or channel maps and 
adds an unnecessary degree of complexity. The signifi- 
cance of the lack of an H I halo has yet to be determined. 
Having such a small number of galaxies observed to this 
depth and modeled using these methods makes it impos- 
sible to extract any reliable trends. This will be remedied 
in the near future via the HALOGAS survey and EVLA 
observations involving some HALOGAS team members. 
For now, we can only compare the results for NGC 4244 
with those of previously modeled galaxies. 

Extended neutral halos have been detected in several 
galaxies listed below. The full resolution la rms noise 
per channel for each external galaxy is given in mJy 
bm -1 . An extended ne utral halo is seen in the Milky 
Way (Levine et aL||2008|), wit h a scale he ight of 1.6 kpc 
dMarasco fc Fraternali||2010[ ) NGC 891 ( |Swaters etaL 



1997 



0.09 mJy bm 1 , 23.4 x 16" resolution) with a scale 
heignt of 1.25-2.5 kpc QOosterloo et al.|2007|), NGC 4559 
with a maximum scale height of 4 kpc ( Barbieri et al. 
2005 1 n " ™ T - 1 ™~ 1 



0.52 mJ y bm 1 , 12.2x24.5" reso lution), as well 
asTTGC 7321 |Matthews fe Wood||2003l ) with a FWHM 
of 3.3 kpc; 0.36-0.40 mJy bm" 1 , 16.2x15.7" resolution 



Uson & Matthews 2003). These are all substantially 
greater than our estimated scale height for NGC 4244 of 
565 pc, although that is only for a thickened disk, rather 
than a disk-halo combination as seen in those galaxies 
listed above. NGC 4559 and UGC 7321 are discussed 
below. 

NGC 4559 has a star formation rate; 0.69 M Q yr" 1 
Heald et al.||2011[ ), approximately 5 times that of NGC 
244. R oughly 10% of the to tal Ft I mass is contained in 
its halo ( jBarbieri et al.|2 005). If H I halos are due to star 
formation, then this would be consistent with NGC 4244 
having a smaller halo, or as found in this work, none at 
all. 

UGC 7321 is also a galaxy with low star formation, 
with a rate of 0.15 M yr _1 . This rate of star formation 
is comparable to that of NGC 4244, and yet a H I halo 
is detected in UGC 7321, but not NGC 4244. At first 
glance, the two galaxies appear quite similar. Both dis- 
play thickened disks with war ps in H I, and by c omparing 
with results for UGC 7321 in |0'Brien et al.| ( |2010j ), their 



rotation curves indicate similar masses, and thus similar 
gravitational potentials. Finally, neither has close com- 
panion galaxies. In spite of these similarities, extended 
wing structures are s een in the vertical profiles shown 
in Figures 2 and 5 of |Matthews fc Wood| ( |2003[ ) , which 
are absent in analogous plots ot 1NLKJ 4244 (not shown), 
drawing a clear and fundamental distinction between the 
morphologies of the two. The presence of a halo in UGC 
7321, a galaxy with a comparable star formation rate 
to that of NGC 4244, indicates that H I halos are not 
necessarily due solely to star formation. 

6.4. Analysis of Distinct features 

It is not the primary intent of this paper to perform 
an exhaustive study of small-scale features. NGC 4244 
is nearly devoid of prominent, energetic H I shells as well 
as expanding features akin to those seen in more actively 
star forming galaxies. However, a small number of no- 
table features are detected (§5) and will be discussed. 

6.4.1. The Shell-Like Feature 

Firstly, we consider the shell-like feature in the ap- 
proaching half. This feature is above and slightly off- 
set radially from a region of star formation (Figure 11 ) 



which extends between 3.5' and 4.6' along the major axis 
in the Ha and 24 /zm data. The feature's proximity to 
a region of star formation as well as the lack of nearby 
external features indicate an internal origin. To estimate 
the energy required to produce it, the number density 
must first be found. For the calculations below we as- 
sume this feature is indeed a shell created by multiple 
supernovae in the disk. Follow ing the approach taken by 
Rand & van der Hulst (1993), we measure the peak flux 
in the limbs ot the presumed shell to acquire a column 
density. To obtain the number density, we then consider 
the path length along the line of sight (I) to be: 



2\/d? + 2dr 



(2) 



Where d is the shell wall thickness and r is the inner 
radius of the shell. Since the shell walls are not well 
resolved, we set d equal to the beam resolution parallel 
to the disk of about 17". From this, a value of 1280 pc 
is obtained for I, which yields a number density of 0.15 
cm~ 3 . 

Now using Figure |12| an estimate for the expansion 
velocity (Kxp) is 25 km s _1 . Together with the estimate 
for the shell's radius in parsecs (i? s heil = 750 pc) and 
assuming a constant expansion, its age is approximately 
2.9 x 10" years. 

We now c alculate the energ y required to produce such 
a shell from Chevalier (1974) where no is the Hi num- 
ber density in cm" 3 , -R s heii is the radius of the shell in 
parsecs, and V cxp is in km s _1 as described above: 



E E = 5.3 x 10 43 n,n 12 RZ 



shell *exp er S 



(3) 



Given these estimates, we find that a single burst of en- 
ergy of roughly 5.4 x 10 53 ergs or the equivalent of 540 
supernovae would be required to produce this shell, (as- 
suming an energy of 10 51 ergs per supernova) If the shell 
instead for med from a continuous sup ply of energy as 
described in McCray & Kafatos ( |1987| ): 

(4) 



E c = 1.16 x 10 4i n o i? s \ cn if J erg 
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where £7 is the age in units of 10 7 years, this would only 
require 1.7 x 10 53 ergs or 170 supernovae. 

Now considering the Ha luminosity of this region we 
ag ain use the total Ha luminosity for NGC 4244 found 
(fl999J). 



Hoopes et al 



The fraction of the total Ha 
emission originating trom this region is approximately 
4%. Thus, its luminosity, uncorrected for extinction is 
found to be 4.0 x 10 38 erg s _1 . 

To correct for extinction, we examine the 24 /im data. 
Summing the flux over the same region as for the Ha data 
and once again assuming a monochromatic flux, the 24 
/im luminosity is approximately 4.8 x 10 39 erg s 



Using the method described in Kennicutt et al. ( 2007 ) , 
we determine the star formation rate, this time with a = 
0.038 and replacing Ltir with the 24 /im luminosity to 
account for localized regions. 

The star formation rate for this region is found to be 
0.0046 M Q j/r _1 , which would account for 3.8% of the 
ongoing star for mation in NGC 4244. However, |Kenni-| 
Icutt et aLll2007l caution that this conversion is intended 
for entire galaxies rather than localized regions, which 
limits its physical meaning. 

Using this star formation rate (assuming it is constant) 
as well as the estimate d age for the shell, a calcu lation 



akin to that presented in McKee & Williams ( 1997[) yields 
an estimated 650-700 supernovae to have gone oft in the 
region, rendering supernovae a very plausible source. 

6.4.2. Curved Feature 

Attending now to the cur ved feature in the reced- 
ing half shown in Figure |13| we examine its potential 
connection to star formation. It is likely that there is 
some correlation between this feature and the underly- 
ing star forming region seen in the both Ha and 24 /im 
data. The best resolved emission is seen in the Ha data, 
which shows a curious ext ensi on of emission from possi- 
bly blown out gas (Figure 13 bottom panel, highlighted 
in boxes) which to some extent fills in the gap defined by 
the H 1 contours, indicating that the H 1 gas may have 
been forced outward in that direction. This structure 
does not appear in the 24 /im data (not shown). When 
examining individual channel maps, it can be noted most 
prominently at 294 km s _1 and 298 km s _1 , a gap in the 
H I forms, once again suggesting a channeling of the gas 
away from the midplane and into higher z in this area 
(Figure [T ) . This hole in the H 1 emission may be analo- 
gous to those seen in other galaxies such as NGC 4559 
(Barbieri et al. 20051. However, it is possible that the 
features are coincident al . 

6.5. Comparison to Previous NGC 4^44 Models 

Our results concerning NGC 4244, using direct visual 
comparisons between models and data, represent an im- 
pr ovement over those obtained using the methods found 



Oiling ( 1996 ) . Firstly, the warp resulting from models 
presented in this paper is more pronounced than that 
previously indicated. Furthermore, the inclination in 
central parts of the galaxy is higher (88 ° as opposed to 
84.5°). The vertical profiles in Figure [5] show a compar- 
ison between the models. It should be noted that in this 
case, the lower inclination in the Oiling model does not 
create the expected rounding near the top of the vertical 
profile and there is a lack of H I in the wings. How- 
ever, this is because the H 1 distribution near the center 



of the Oiling model is much higher than in the models 
presented here. This is also seen in bv diagrams of the 
Oiling model (not shown) . Additionally, there is no need 
for a flaring layer in the current models in order to fit the 
data, although one may be present. These discrepancies 
cannot be accounted for by improved observations alone 
and appear to result from differences in technique. 

7. SUMMARY AND CONCLUSIONS 

We present tilted ring models based on deep 21-cm 
observations of NGC 4244. In these observations, we 
note a thickened H 1 disk, with thinning at outer radii 
and evidence for a warp. We also note asymmetries in 
the approaching and receding halves as well as localized 
features. The models tested, in order of approximate 
increasing complexity, include a single disk with constant 
inclination, a single disk with a warp component along 
the line of sight, a single flaring disk, a thin and thick 
disk, and combinations of these. 

From these models, we conclude that the disk of NGC 
4244 is substantially warped both perpendicular to as 
well as along the line of sight. We do not detect an 
extended Hi halo, but with a scale height of 0.6 kpc, 
the disk is quite thick. If H 1 halos are due to disk-halo 
flows, then given the low star formation rate in NGC 
4244 compared to other galaxies with H 1 halos, the lack 
of a halo is unsurprising. However, it is still to early to 
establish a trend, which will hopefully become apparent 
when additional cases are considered. 

In spite of not detecting an extended halo, we do detect 
a vertical lag in rotation speed of —9^?, km s _1 kpc -1 
and — 9±2 km s _1 kpc -1 in the approaching and receding 
halves respectively. The magnitude of this lag decreases 
outward with radius to — 5± 1 km s _1 kpc -1 (approach- 
ing) and —4 ± 2 km s _1 kpc -1 (receding) near 10 kpc, 
although subtle changes are within the error estimates. 
Beyond 13 kpc, the S/N ratio no longer allows for a reli- 
able lag estimate, so we cannot discern whether the lag 
vanishes entirely at some point. 

This is potentially the first instance where models favor 
a lag in Hi in a galaxy where an H 1 halo does not exist . 
There may be a lag in UGC 1281 ( |Kamphuis et al.|2011[ ), 
in which a halo is also absent, but in that case the lag 
cannot be distinguished from a warp component along 
the line of sight without deeper observations. Given fluc- 
tuations in the position angle seen in the zeroth and first 
moment maps, it is possible that the warp in NGC 4244 
begins closer to the center than in the models presented, 
which may indeed lessen or eliminate the need for a lag. 
However, such a scenario is unlikely as these fluctuations 
may instead be due to localized features, a bar or spiral 
arms, which would be more in-line with current under- 



standing of warps (Briggs|1990||van der Kruit & Freeman 
20111 



Modeling aside, we observe two substantial localized 
features which may correlate with star formation: a shell- 
like feature and an unexplained curved structure origi- 
nating in the midplane and possibly connected with a 
feature extending up to a height of 1.5 kpc. A third dis- 
tinct feature, a faint H 1 streak is anomalous. Finally, no 
additional substantial energetic features are detected, al- 
though numerous small-scale shell-like features are seen, 
but not displayed. 
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We hope to use information gathered from our models 
of NGC 4244 in addition to information from the rest 
of the HALOGAS galaxies to further the development 
of trends concerning H I kinematics, warps, extra-planar 
gas and any connections with disk-halo interactions, as 
well as accretion. 
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